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Abstract. We prove dispersive decay estimates for the one-dimensional Dirac operator and 
use them to prove asymptotic stability of small gap solitons in the nonlinear Dirac equations 
with quintic and higher-order nonlinear terms. 

1. Introduction 

Asymptotic stability of solitary waves in the nonlinear Schrodinger equation has been consid- 
ered in the space of three dimensions with a number of analytical techniques [20|, 127 } [2^ 1 [8j \6[ rH] . 
Only recently, the asymptotic stability of solitary waves was extended to the space of two di- 
mensions [HI Q2] and one dimension O Q2] . 

Relatively little is known about the asymptotic stability of solitary waves in the nonlinear 
Dirac equations, which can be considered as a relativistic version of the nonlinear Schrodinger 
equation. Asymptotic stability of small bound states in the nonlinear Dirac equations in three 
dimensions was constructed by Boussaid [2j. Global existence and scattering to zero for small 
initial data were obtained by Machihara et al. \15\ I16j. also in the space of three dimensions. 

We shall consider the asymptotic stability of solitary waves in the nonlinear Dirac equations 
in one dimension. Since the energy functional of the Dirac equations is sign-indefinite at the 
linear wave spectrum, it is generally believed that the solitary waves (referred to as gap solitons) 
must be energetically (and nonlinearly) unstable. Indeed, gap solitons are more disposed to 
spectral instabilities in the sense that unstable eigenvalues may exist in a large subset of the 
existence domain [4] . However, the limit of small gap solitons corresponds to the nonrelativistic 
limit, when the nonlinear Dirac equations can be reduced to the nonlinear Schrodinger equation 
|16j . In this limit, when the cubic nonlinear terms are considered, the gap solitons in one 
dimension are typically stable both spectrally and orbitally. It is hence an interesting question 
to study the nonlinear asymptotic stability of the spectrally stable small gap solitons. 

The spectral information is difficult in the case of the homogeneous nonlinear Dirac equation 
even in the limit of small gap solitons. Isolated nonzero eigenvalues and resonances at the 
end points of the continuous spectrum occur commonly in the problem [US]- To simplify the 
spectral information, we add a bounded exponentially decaying potential to the one-dimensional 
nonlinear Dirac equations and consider a local bifurcation of the small gap solitons from an 
isolated eigenvalue of the self-adjoint Dirac operator. In this way, our approach is similar to 
the one used by Mizumachi [T7] for the nonlinear Schrodinger equation and by us [12] for the 
discrete nonlinear Schrodinger equation (see also [7] for similar results). 

We shall avoid the dispersive decay estimates in weighted I? spaces, which are difficult for 
the nonlinear Dirac equations (in contrast with the nonlinear Schrodinger equations). We shall 
instead derive the Strichartz estimates directly from the Mizumachi estimates. The balance 
between Strichartz and Mizumachi estimates allows us to control both the nonlinear terms and 
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the modulation equations for small gap solitons and thus to prove their asymptotic stability 
for the nonlinear Dirac equations with quintic and higher-order nonlinear terms. 

The article is organized as follows. Section 2 introduces the nonlinear Dirac equations. 
Section 3 contains information about the small gap solitons. Section 4 reports on linearization 
and spectral stability for small gap solitons. Section 5 derives the modulation equations for 
parameters of gap solitons as well as the time evolution equation for the dispersive remainder 
term. Section 6 describes the spectral theory for the one-dimensional Dirac operator. Section 7 
deals with the linear dispersive estimates for the semi-group associated with the Dirac operator. 
Section 8 gives the proof of the main theorem. 

We finish this section with the list of useful notations. 

The inner product for complex- valued functions in L 2 (M) is denoted by 



(1) 

For any / G L 2 
(2) f(k) ee F{f) : 



(f,9)^ ■= / f(x)g(x)dx. 



, we define the Fourier transform and its inverse by 

-i poo i roc 

1 ' " - -ixk - - 1 - - ' 



f(x)e-™"dx, f(x) = T-\f):-- 



J- 



f{k)e ixk dk. 



Sobolev spaces are denoted by W S < P (R) for s > and 1 < p < oo so that H S (R) = W S < 2 (R) 
and L P (R) = W° ,P (M). Beside Sobolev spaces, we will use Strichartz spaces Jj^lJx cind L%L^ 
defined for 1 < p, q < oo by the norms 



(3) 



i/p 



ljf(x,-)\\i Pt d X y q , 



where T > is an arbitrary time including T = oo. 

Notation (x) = (1 + x 2 ) 1 / 2 is used for the weights in L% norms. The constant C > is a 
generic constant, which may change from one line to another line. A ball of radius 5 > in 
function space X centered at £ X is denoted by B$(X). 

Pauli matrices are defined by 



The 2-by-2 identity matrix is denoted by Id. 

Scalar functions are denoted by plain letters and vector functions with two components are 
denoted by bold letters. For clarity of notations, we do not write the second arguments for 
W S,P {U), H S {E), and L 2 (M) when it is used for scalar or vector functions. 
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2. The nonlinear Dirac equations 
Consider the nonlinear Dirac equations 

i(u t + u x ) + v = duW(u, v), 



where (x,t) G M 2 , (u,v) G C 2 , and W(u,v) : C 2 - 
following three conditions: 

• symmetry W(u, v) = W(v,u); 

• gauge invariance W(e ie u,e te v) = W(u, v) for any 9 6 

• polynomial in (n, v) and (n, v). 



i(vt - v x ) + u = dvW(u, v), 

R is a nonlinear function which satisfies the 
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A general expansion of the nonlinear function W(u, v) satisfying the three properties above 
starts with quadratic and quartic terms 

(5) W = (3(x)(\u\ 2 + H 2 ) + j(x)(uv + uv) + W N (u, v), 
where /3(x), r y(x) : M — > M are bounded and decaying potentials, 

(6) Wn = oi(|ii| 4 + \v | 4 ) + a2|n| 2 |t>| 2 + a^(uv + uv) 2 + a4(|u| 2 + |v| 2 )(-ut> + uv) 

is the nonlinear (quartic) potential, and (01,02,03,04) G M are numerical coefficients. 

The standard example of the nonlinear term occurs in the context of Bragg gratings, where 
P(x) and 7(x) model optical defects in the periodic grating, whereas 

(7) W N = o(M 4 + 4|u| 2 H 2 + H 4 ), ael 

models the nonlinear coupling terms \10\ [TT] . 

Another example is relevant to the massive Gross-Neveu model for spinors in relativity 
theory [I], 



(8) 



Wn = a(uv + uv) x 



a G 



In other applications, Wn may start with terms of the sixth and higher orders. The follow- 
ing nonlinear potential is derived in the context of the Feshbach resonance for Bose-Einstein 
condensates [2JJ, 



(9) W N = a(\u\' z + \v\' z )\u\ 2 \v\ 2 , oG 

Let us introduce the 2-by-2 Dirac operator in one dimension 

-id x + P(x) <y(x)-l 



(10) 
where 
(11) 



D 



j(x) 

-id x 
-1 



1 



-1 

id x 



id x + fi{x) 



V{x) 



D + V(x), 



P(x) j(x) 
j(x) p{x) 



The nonlinear Dirac equations can be rewritten in the abstract evolutionary form 



(12) 



. du 

'~dt 



^u + N(u), N(u) = V„Wn(u,v), u 



u 

V 



Vq 



where N(u) = C(||u|| 3 ) as ||u|| — > in any norm that forms Banach algebra (e.g. in H S (R) for 
s > 4). For the potentials ([7]) and (jHJ), we have explicitly 



N(u) = 2o 



(\u\ 2 + 2\v\ 2 )u 
(2\u\ 2 + \v\ 2 )v 



N(u) = 2o 



\v\ 2 u + v 2 u 
\u\ 2 v + u 2 v 



For the technical reasons, these cubic nonlinear functions are not sufficiently small when (u, v) 
decays to zero. As a result, we shall consider a more general class of the homogeneous polyno- 
mials of Wn (of even degree). Our arguments will be valid for the quintic nonlinear functions 
which are generated from the polynomial Wn of degree six, e.g. from the function ([9]). 
Local existence of solutions in Sobolev space can be proved with standard methods |10j . 



Proposition 1. Let uq G H S (WL) for a fixed s > \ and assume that W satisfies the three 
conditions above. There exists a T > such that the nonlinear Dirac equations admits a 
unique solution 

u(t) G C([0,T],H s (M))nC\[0,T},H s -' 
where u(t) depends continuously on the initial data u(0) = Uq. 
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If the nonlinear functions Wn(u,v) depends on \u\ 2 and \v\ 2 only, e.g. for (J7J) and Q, global 
well-posedness in H S (M) with s G N can be proved [10]. Little is known about the global 
solutions even for small initial data for the general nonlinear Dirac equations, e.g. for ([8]). 

3. Stationary small gap solitons 

Under the assumptions that f3(x),j(x) G L°°(IR), Dirac operator % is a densely defined, 
self-adjoint operator in L 2 (M) with the domain H 1 ^). We shall further assume that 

/3(x), j(x) — > as |x| — > oo 

at an exponential rate. The potentials represent a relatively compact perturbation to the 
unbounded differential operator. By Weyl's Theorem, the spectrum &{%) C M contains the 
continuous spectrum at 

ff c (H) = (-oo,-l]U[l,oo). 

To simplify the construction of stationary small gap solitons, we assume that H admits only 
one simple isolated eigenvalue in the gap (—1,1) of the continuous spectrum a c {H) and no 
resonances at ±1. Hence we add the assumption. 

Assumption 1. Assume that 

• /3,7 G L°°(1R) and there is C > and k > such that 

\P(x)\ + \-y(x)\ < Ce~ K ^, iel. 

• a(rl)\a c ( / H) = {ujq}, where ujq G (—1,1) is a simple eigenvalue of H with the L 2 - 



normalized eigenfunction uq 6 H 1 ( 
• No resonances occur at the end points ±1 of cr c (1-L). 

Stationary gap solitons are given by 

(13) u(x,t) = U(x)e- iLut , v(x,t) = V{x)e- iu,t , 

where w S f is a parameter and U = [U, V] T G C 2 satisfies the system of differential equations 

(14) (H - cjI)V + N(U) = 0. 

If U G i? 1 (R), then U G C(R) and U(x) as \x\ -)■ oo thanks to Sobolev's embedding 
of fl' 1 (R) to C$(R). By Lemma 3.1 in [4], the stationary solution U G fl' 1 (R) satisfies the 
symmetry 

(15) U(x) = V(x), x£R. 

For the example of nonlinear function ([7]) with a = when no linear potentials are present 
f5{x)^{x) = 0, the stationary gap solitons are given in the explicit form for any cj G (—1, 1) 

(16) U(x) = V ^ = 

\/l — cj cosh(y 1 — oj 2 x) + + wsinh(v 1 — uj 2 x) 

In particular, ||U||l°o — > as uj — > —1, which indicates the limit of small gap solitons. 

As we explained in the introduction, the spectral information is difficult in the case of 
homogeneous Dirac equations (without decaying potentials). If /3(x) and 7(x) are nonzero 
and Assumption Q] is used, the stationary gap solitons are not known in the explicit form but 
the local bifurcation technique allows us to find a family of small gap solitons in a one-sided 
neighborhood of to = ujq. To make it more precise, let us assume that the nonlinear function is 
a homogeneous polynomial in its variables. 
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Assumption 2. Assume that 

N(aU) = a 2p+1 N(U), a G R, 

for a fixed integer p > 1 . 

Proposition 2. Lei Assumptions [7] anaTJI &e £rue and 

(17) (u ,N(u )) L 2 >0. 

For sufficiently small e > 0, £/iere is a family of solutions U G of system Jj^p /or any 

u> G (wo,w + e) suc/i f/iai i/ie map (u>o,w + e) B uj h- >• U G £T (R) is C 1 and 

(18) ||U-ou ||ifi = 0(a 2p+1 ), |w - wo | = 0(a 2p ), as a -> 0. 



Proof. Thanks to Assumption [TJ we use the decomposition 

U = au + V, a G R, (u ,V) L 2=0. 

Let Po : L 2 (R) — > Ran(% — ojqI) C L 2 (R) be the orthogonal projection operator so that 
V = PqV G Ran(% — ujqI). The stationary equation (|14p becomes the following system of two 
equations 

f P (H-uI)P V + P o N(au + V) = O, 
\ (uj -oo)a+ (u ,N(au + V)) £2 = 0. 

Operator P Q (H - uj I)P : F^R) -> L 2 (R) is invertible for w near wq, whereas N(U) is a 
(j2p+i f unc tion near G i? 1 (M). By the Implicit Function Theorem, there is a unique C 2p+1 
map R 9 a i— )• V G -fT 1 (R)nRan(Lo — w o-^) such that V satisfies the first equation of the system 
and there are a > and C > such that for all a G (— Oq, a ), 

(19) ||V - a 2 P +1 P (?{ - woir^oNK)!!^ <Ca 4 ? +1 . 

Let us substitute the map KBa^VG F 1 (R) n Ran(Lo — u)qI) to the scalar equation 

F(a,w) = (w -w) + a" 1 (uo,N(au + V)) i 2 =0. 

Thanks to the bound (|19p and Assumption [2j there are oq > and C > such that for all 
a G (— Ooj a o); there is only one solution of F(a,oj) = for uj = uj{o) satisfying the bound 

(20) \uj + a 2p (u , N(u )) L 2 - u\ < Ca Ap . 

Under condition (fTTl) . we have u > loq and the bounds (fl8l) follow from (fl9j) and (|20l) . □ 
Remark 1. Proposition^ is valid if 

(21) (u ,N(u )) L 2 <0, 

ou£ the family of solutions U G H l {E) of system (TT^P exist for oj G (wo — e,w) under the 
condition [2l\). 
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4. Linearization and spectral stability 
Linearization is performed after writing 



(22) 



u(x, t) 
v(x, t) 



-iwt 



-iwt 



U(x) + U 1 (x)e xt + U 2 (x)e xt 
V(x) + V 1 {x)e xt + V 2 (x)e xt 



and neglecting quadratic terms with respect to the vectors 

Ui = [U u Vxf G C 2 , U 2 = [U 2 , V 2 f G C 2 . 

From the derivatives of N(U), we obtain the expansion in any norm that forms a Banach 
algebra with respect to the pointwise multiplication, 

(23) N(U + Ui) = N(u) + y 11 u 1 + y 12 u 1 + o(||u 1 || 2 ), 

where Vn and V± 2 are 2x2 matrices with exponentially decaying coefficients, which are given 
explicitly by 



(24) 



ii 



&vu w » d lv w » 

4u W N 4v W N 



11) 



V5 



12 



sfow N d^ v w N 

4u W N 4v W N 



12- 



Substituting (|22j) to the nonlinear Dirac equations (jl2j) and using expansion (|23j) . we obtain 
the linear eigenvalue problem 



(25) 



uui = {%- w/)Ui + FnUi + vi 2 u 2 , 

-iAU 2 = (Ho - w/)U 2 + V12U1 + V n U 2 . 



We should distinguish the self-adjoint operator H u : H l { 



L 2 (JH) given by 



H 



n-ui 

Hq-ljI 



V n Via 
Via Vn 



and the non-self-adjoint linearization operator L, 
(26) 



^ L ! 

-iaH, : tf 1 



L 2 (R), where 



a 



I 
-/ 



Both operators act on [Ui,U 2 ] T . 
Symmetry (|15p imply that 



d 2 Du w N = d 2 Vv w N , d 2 m w N = d 2 vv w N , d 2 uv w N = d 2 DV w N . 

By Theorem 4.1 in [4], the self-adjoint operator H u and the linearized operator L w can be 
block-diagonalized. Let S be an orthogonal matrix given by 



S 

Direct computations show that 

(27) S-t-HuS 

(28) S^aH^S 





/ 1 





1 





\ 


1 





1 





-1 




72 





1 





1 






V 1 





-1 





/ 



H+ 


/ 
—I 








H- ' 
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where H± are two-by-two Dirac operators given by 

-id x + (3(x) -oo ±( r r(x) 



(29) 



1) 



+ V ± (x), 



±( 7 (s) - 1) id x + P{x) - u 
and V±(x) are 2-by-2 matrices with exponentially decaying coefficients given by 



(30) 



V+ 



d 2 Qu w N ± d 2 QV w N d 2 QD w N ± d 2 Qv w N 

d 2 JU W N ± ^Wjv ± 



Thanks to the symmetry of the nonlinear Dirac equations @ with respect to the gauge 
transformation, the linearized operator L u has a nontrivial kernel because 

U 



(31) F = t ~j G Ker(L w ) = Ker(H u ), 
or explicitly 

(32) (H - uI)XJ + V n TJ - V 12 V = 0. 
The eigenvector F generates a two-dimensional generalized kernel 

U 



span{F,G} C N g (L u 



-a 



U 



(33) 

such that L^G = F, or explicitly 

(34) (H - ojI)8 u U + Vud u U + V X2 dJJ = U. 

The Jordan block is two-dimensional, that is, no H g fl^QR) solving L W H = G exists, if 



(35) 



dco ^ 



L 2 



^0. 



Constraint (|35p is satisfied for small a in Proposition [2] under condition (|17p . In the same limit, 
the spectra of the linearized operator L w and the self-adjoint operator H w are characterized in 
the following proposition. 

Proposition 3. Let assumptions of Proposition^ be satisfied. For sufficiently small e > and 
for any to G {ojq,uiq + e), we /jaue 

o-(iL w ) = (-oo, -1 - w] U (-oo, -1 + w] U {0} U [1 - oo, oo) U [1 + w, oo) 

and 

cr(iT + ) = (-oo,-l -w]U{a;i}U [l-w,oo), <r(.ff_) = (-oo, -1 - u] U {0} U [l-u,oo), 

where oo% = O(\oo — too\). The zero eigenvalue is double for L w and simple for whereas the 
eigenvalue to\ is simple. No resonances exist at the end points of the continuous spectrum of 
L w and H± . 

Proof. The proof holds by perturbation theory. The self-adjoint operator H + is represented by 

H + = U-uoI + V + , 

where ||V+||£,°° = O{o? p ) and \co — ooq\ = 0{a 2p ) as a — > (parameter a is used in Proposition 
[2]). By Assumption [TJ no resonances exist in T~L and, therefore, no new eigenvalues bifurcate to 
the gap (— 1 — to, 1 — to) of the continuous spectrum of H + from the non-resonant points ±1 — to 
for small a > 0. By the perturbation theory, the only eigenvalue at for a = becomes the 
eigenvalue to\ = O(o? p ). 

The self-adjoint operator iJ_ is given by 

H- = a 3 Hcr 3 -ojI + 
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where ||VL||ioo = 0(a 2p ) as a — > 0. The same perturbation theory applies to self-adjoint 
operator except of the fact that the only eigenvalue at for a = is preserved at for 
a > thanks to the gauge invariance, which results in the exact relation 



U 

-V 



0. 



Similarly, the double zero eigenvalue of L u is preserved at by the gauge invariance as the 
generalized kernel (133H . whereas the continuous spectrum does not lead to resonances at the 
end points or to new eigenvalues for small a > 0. □ 

5. Projections and modulation equations 
By Proposition [3j we have 

(36) NgiL^) = span{F,G}. 
Recalling matrix a from (|26p . we obtain the adjoint operator 

L* = (-iaH w f = iH*a* = iH u a, 
which has the generalized kernel 

(37) N g (Ll) = span{<rF,aG}. 

Any vector [Ui, U 2 ] T in the invariant subspace of the linearized operator in L 2 (R), which 
is an orthogonal complement of the generalized null space N g (L w ), has to satisfy the symplectic 
orthogonality conditions 

loo) / (U,U 1 ) i2 + (U,U 2 ) i2 =0 ) 

{ 1 1 (djLi,Ui) L 2 - (d u \j,v 2 ) L 2 = 0. 

If U2 = Ui, the symplectic orthogonality conditions (I38p can be rewritten in the explicit form 

(39) Re(U, V X ) L 2 = 0, Im^U, Ui) L 2 = 0. 

Using symplectic orthogonality conditions, we now set up modulation equations for nonlinear 
dynamics of small gap solitons. By Proposition [H we have at least local solutions of the 
nonlinear Dirac equations (|12p . Now we look for local solutions in the form 

u(x, t) = e- i9 V [[/(x;o;(t)) + [/i(x,t)], 
v(x,t) = e" ie W [V(x;u{t)) + Ui(x,t)] , 

where we write explicitly the dependence of the stationary solution U = [U, V] T on uj. The 
time evolution problem for Ui = \U\, Vi] T is given by 

(41) i^i = (U- uI)Vx - iuduU - - w)(U + Ui) + N(U + Ui) - N(U). 

Using the symplectic orthogonality condition (139)) on U, we obtain the modulation equations 
on u(t) and 9{t): 



(40) 



(42) 
where 



wRe(S w U,U-Ui)£a + (0 -w)Im(U,Ui) L 2 = F u 
uIm(d 2 V, Ui)^ + - w)Re(9 U) U, U + Ui) L 2 = F 2 , 



Fx = Im[(U,N(U + U 1 )-N(U)) i2 + (U 12 U-U 11 U,U 1 ) L2 ] , 

F 2 = Re[(a w U,N(U + U 1 )-N(U)) i2 -(U 12 a w U + Una tJ U,U 1 ) i2 ] 
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and equations (|32p and (|34p have been used. The following result shows that the right-hand 
side of system (l42j) is quadratic with respect to the perturbation vector Uj. 

Proposition 4. Let assumptions of Proposition be satisfied. Fix small e > and 5 > For 

any u G (wo,wo + e) and any Ui G Bs(L°°), there is C €) $ > suca iaai 

(43) l^il + l^l < C e , 5 <U2,U> i2 . 



Proof. We use the expansion (I23p . constraints on matrices (|24p and elementary properties of 
inner product ([I]) to show that the linear terms in Ui vanish in the expression for F\ and F^. 
For instance, F\ has the following linear terms in Ui: 

Im [(U, VnUi + Fi 2 Ui) L 2 + (VisV - FuU,Ui) L2 ] 
= Im [(^U,U!) L2 - {V n V,Vi)v + <v£u,Ui) L 2 + <yi 2 U,Ui) L2 ] = 0. 

Similar computations holds for linear terms of F 2 . Together with smoothness in Assumption 
[2j this computation shows that both terms Fi, F 2 are quadratic in Ui in the sense of (|43p . □ 

Setting Ui = Ye ,e , we rewrite the time-evolution equation (I4ip in the equivalent form 
dY 

(44) i— =HY + e- ie F, F = -iwd w U-(0-a;)U + N(U + Ye tfl )-N(U). 
We are now ready to formulate the main theorem of this article. 

Theorem 1. Assume Assumption^ Assumption^ with p > 2, and condition fli?| ). Fix e > 
and <5 > sufficiently small such that 8(0) = 0, w(0) G (wo,wo + e); and Y(0) G B$(H l ). There 
exist eo > e, 6*00 G R, Woo S (wo. w o + ^o), (oj,8) G C 1 (M + ,1R 2 ), and 

Y(i) G C(R + ,H 1 ) n L 4 (M + , F°°) 

suca that (u),6)(t) solve the modulation equations (4*fy > Y(t) solves the evolution equation \44\ )> 
and 

lim I — / uj(s)ds I = ^oo, lim = oJqo, lim ||Y(t)||x°o = 0. 

t— >oo \ Jg y t— >oo t— >oo 

We shall prove this theorem in the remainder of the article. To do so, we shall develop first 
the spectral theory for the Dirac operator % and obtain the dispersive decay estimates for the 
semi-group e~ ltH acting on the continuous spectrum of % . 

6. Spectral theory for operator H 

Let us consider the spectral problem rlu = Au or explicitly, 

/ 45 ) | -iu'(x) + (3{x)u{x) + (7(1) - l)v(x) = Xu(x), x G R> 

\ iv'(x) + f3(x)v(x) + (7(2;) — l)n(x) = Xv(x), 

Recall that 

ff c (H) = (-oo,-l]U[l,oo). 

Here we develop the scattering theory of wave operators for the Dirac operator T~L. A similar 
theory for the Schrodinger operators on an infinite line goes back to the works of Weder |25tl26| 
and Goldberg & Schlag [9]. 
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Let us first define the Jost functions for A G (—00, —1] at one branch of a c {%). To do so, let 



us parameterize (—00, —1] by A = — yl + k 2 for k G M and consider solutions of system ([35 
according to the boundary conditions 



(46) 



u (x; k) 



1 

a±(k) 



e ±ikx as x 



±00, 



where a±(k) := k 2 ±k. The following proposition gives the construction of Jost functions. 

Proposition 5. For any there exist unique Jost functions u =t (x; k) such that 



lim 

x— >±oo 



u ± (-;fe)-[l,a+] J e 



T+ikx 



0. 



Moreover, 

• Ifk^O, then u ± (-; k) G L°°(R). 

• If k = 0, then u =t (x; 0) may grow at most linearly in x as x — > =Foo. 

• As k — > ±00, both u + (x;k)i and u~(x;k)i are bounded, u =t (x;/c)2 grows linearly in k, 
and u =1= (x; k)2 decays inverse linearly in k. 

Proof. Setting u^(x;k) = m =l: (x; k)e ±lkx and using the Green function technique, we obtain 
an integral equation for the Jost functions m =l: (x; k) 

-±oo 



(47) 
where 
(48) 



m (x;k) 



1 

a± 



+ 



G (x - y; k)V(y)m (y;k)dy, 



a=p — a±e 



^2ikx 



1 



D ^f2ikx 



1 



a ^2ikx 



3 r 2ikx 



a± — a T e n 

Under the assumption of fast decay of V(x) to as \x\ — > 00, the standard theory gives 
solutions m =t (-; k) G L°°(IR) of the integral equations (|47p for k 7^ and the scattering relation 
between the solutions for all x G M including x —> =Foo 

m + (x; k) = a + (k)m~ (x; —k) + b + (k)mT (x; k)e~ 2%kx , 
m~(x; k) = a~ (k)m + (x; —k) + b~ (k)m + (x; k)e 2lkx , 



(49) 
where 



a ± (k) 
b ± (k) 



1 ± 



2ik 



[a T [V(x)m ± (x;k)]i + [V r (x)m ± (x; k)} 2 ) dx, 



2ik 



a±[V(x)m (x;k)]i + [V(x)m (x;k)]2) e 



+2ikx 



dx. 



It follows from the explicit expressions for a^{k) and b ± (k) that 



(50) 
where 
(51) 



a ± (k) 



± 



2ik 



, b ± (k) 



=F7 



as A; — > 0, 



7 ± = / (/3(x) + 7(x))(m^(x;0) + m^(x; 0))dx. 
If = 0, the Jost functions m ± (x;0) satisfy the integral equation 



(52) 



m ± (x;0) 



" 1 " 




1 


+/ 

J X 



/>±oo 

/ G ± (x-y;0)y(y)m ± (y;0)d 2/ , 
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where 

-fcfc, 



(53) G (x;0) 



X + l X 

x x — i 



Thanks to the fast decay of V(x), existence of locally bounded function m^(x; 0) follows again 
from the standard theory. The linear growth of m ± (x; 0) as x — > =Foo follows from the integral 
equations (f52|) if 7 =t ^ 0. 

Finally, as k — > ±oo, a± grows linearly in k, decays inverse linearly in k, whereas G ± (x; k) 
remains bounded. The asymptotic behavior of solutions m ± (a;; k) of the integral equations (|47p 
follows the asymptotic behavior of the limiting functions (|46p in k as k — > +oo. □ 



Remark 2. Proposition^ eliminates the possibility of embedded eigenvalues in the continuous 
spectrum a c (H) because the space of solutions of the Dirac system |^5|) for A < — 1 is spanned 
by the two fundamental solutions u^(x;k) with no decay to zero as x — )• ±oo. 

The following proposition summarizes the relations on the scattering coefficients in the scat- 
tering relation fj49f) . 



Proposition 6. For any k £l, we have 

(54) a + (k)=a~{k), b + (k) = -b~ (-k), 



,/TTp — h_ 

(55) a+(-k) = a+(k), b + (-k) = -y==——b + (k), 

V 1 + k + k 

and 

(56) |a+(£;)| 2 = l + ^E|^|&+(A0^ 

V 1 + k z + k 



(57) 



Proof. Inverting the scattering relation (|49p . we obtain the constraint on the scattering coeffi- 
cients for all k G M 

a+(k)a-(-k) + b+{k)b~{k) = 1, 
a+(k)b~(-k) + b + {k)a-{k) = 0. 

Let H /r (ui,U2) denote the Wronskian determinant of any two solutions ui and U2 of the 
Dirac system (|45p . It is clear that iy(ui,U2) is constant in x G R. Therefore, W(ui,U2) can 
be computed in the limits x — > ±00. Using boundary values (j4U|) and scattering relations (|49|) . 
we obtain 

(58) W(u + , u - ) = uf(x; k)u^ (x; k) — u^x; k)u^ (x; k) = —2ka + (k) = —2ka~(k). 

This result together with the second equation of system (I57p gives relations (I54p . The first 
equation of system (I57p implies now for all k 6 M. that 

(59) a + {k)a + {-k) - b + (k)b + {-k) = 1. 



Let (uk,Vk) denote any solution of the Dirac system (j45l) for A = —y/l + k 2 . It is checked 
directly that 

\uk\ 2 -\vk\ 2 and u^ k u k - V-kVk 



are constant in x. Using boundary values ([46 p and scattering relations (|49|) again, we obtain 

(60) (Vl + k 2 + k)(l - \a + (k)\ 2 ) + (Vl + k 2 - k)\b + (k)\ 2 = 0, 

(61) {\fl + k 2 + k)l~{k) + (\/l + k 2 — k)b + {k) = 0. 

These identities together with equation (f^|) give relations (f^5]) and (f5T)j) , □ 
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Remark 3. Identity A56\) shows that \a + (k)\ > 1 for all k G R. This excludes embedded 
resonant states with a + (k) = 0. There is still a possibility of end-point resonances at k = 0, 
since a + (k) is generally singular as k — > 0. We say that the end points ±1 are resonances if 
there exist a solution u G L°°(R) of the spectral problem f or A = ±1. If this is the case, 
then 7 + = and lim^o a + {k) exists. 

We shall now define the Jost functions for A G [1, oo) at the other branch of a c (T~L). Similarly 
to the analysis for A G (— oo, —1], we can parameterize [1, oo) by A = Vl + k 2 for k G R and 
consider solutions of system (j45]) according to the boundary conditions 



v =t (x; k) 



-a± 
1 



e ±lkx as x — > ±oo. 



Using a similar Green's function formulation, Proposition [5] can be extended to functions 
v ± (x; k). In what follows, we will not treat functions v ± (x;/c) for A G [1, oo) but will only 
be working with functions u =t (x;/c) for A G (— oo, — 1]. This approach does not limit any 
generality. Moreover, we note the particularly remarkable case. 

Remark 4. If f3(x) = 0, the Jost functions are related by 



v (x; k) = 

thanks to the symmetry of the Dirac system 



-1 

1 



Let Ry.(X) = {T~L — A/) -1 be the resolvent operator, defined as a bounded operator from 
L 2 (R) to L 2 (R) for any A ^ a{J-L). Using the Jost functions and the standard limiting absorption 
principle, the resolvent operator is extended to the continuous spectrum as a uniformly bounded 
operator from L 2 (R) to L 2 _ Q (R) for any a > ^. Let us denote the limiting operators by 

R±(\) := limR n (X±ie), \£a c (H), 

depending on whether A — > cr c (T-L) from the upper or lower halves of the complex plane of A. 

The following proposition allows us to express -R^(A) for A G a c {T-i) in terms of the Jost 
functions. According to the previous remarks, it is sufficient to consider A G (— oo,— 1]. The 
arguments for A G [l,oo) can be developed similarly. 

Proposition 7. For any A G (— oo, — 1] and any fixed a > i, operators R^(X) ■ -^q(R) | — > 
L 2 _ a (R) can be represented by the integral kernel in the form 

(62) rp±mif T ^- ±:L { u + (x;±k)[o-iu-{y;±k)] T , x > y, 
[b2) [thiWlWV) ~ 2ika+(±k) \ u-(x;±k)[a 1 u+(y;±k)] T , x<y, 

where k < and A = — \f\ + k 2 . 

Proof. Let us consider the solutions of the linear system for a fixed y G R and A = — yl + k 2 , 

(63) (H-\I)[R+(\)](x,y)=6(x-y)Id, 
which satisfy the asymptotic behavior, 

(64) [R+(X)](x,y) r^e**!*-"!, as \x - y\ ^ oo. 

The function [i?^(A)](x, y) decays exponentially as \x — y\ — > oo if k is extended off the real axis 
with Im(fc) > 0. Since Re(A)Im(A) = Re(/c)Im(/c) and Re(A) < —1, we understand that the 
behavior (|64p recovers the limiting resolvent operator i?^(A) defined for Im(A) > if Re(/c) < 0. 
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(65) 



u 


-{ 


V 





For the first column vector of the linear system (j63|h denoted by (u,v), we obtain 

c(y,k)u + (x;k), x>y, 
d(y,k)u~(x;k), x<y, 

where the behavior (|64p is satisfied thanks to the boundary conditions ()46j) . Parameters (c, d) 
are to be determined. 

Matching conditions across the point x = y sets up the linear system for c and d with the 
unique solution, 

C (^) = ^¥% d(y,k)- lU+ ^ 



W(u+,u-)' vy ' ; W(u+,u-)' 

where W(u + ,u~) = -2ka+(k) by identity (158]) . 

Similarly for the second column vector of the linear system f)63[) , we obtain the same expres- 
sion (|65p with a different solution of the linear system for (c, d) , 

c (y, «) = ttt- / _i_ — zr> fc) - 



W(u+,u-)' v y W(u+,u-) 

Using the Pauli matrix a\, we arrive to the expression (|62p for R^(X). The expression for R^(X) 
is found by the replacement of A: by —k. The exponential decay as \x — y\ — > oo occurs now for 
Im(fc) < 0. The limiting resolvent operator i?^(A) is defined for Im(A) < if Re(/c) < 0. □ 

The following proposition describes A-uniform bounds on the limiting resolvent operators 
i?^(A) in weighted spaces. In order to exclude problems at the end points A = ±1, we assume 
that no end-point resonances occur at k = (Assumption [1]). Thanks to Remark [3l it is 
equivalent to assume that 7 + ^ 0. 

Proposition 8. Let 7 + ^ in \51\) . For any a > |, there exists constant C a > such that 

(66) sup ||^(A)|| L ^ i2 <C a . 

\\\>l 

In addition, for any a > 1, there exists constant C a > such that 

4( 



(67) sup ||1J±(A)|| C7 a . 



|A|>1 

Proof. We recall that |a + (/c)| > 1 (Remark [3]). Thanks to the asymptotic expansion (|50p . if 
7 + 0, then ka + {k) ^ for any k € R. Using this result and Proposition [71 we construct 



[K£(\)}(x,y) 



(l + a J)o/2(l + y2)a/2- 

By Proposition ir^-jfc) e L°°(R) for every 7^ and u (x;0) grow at most linearly in 
x as x — > =Foo. Therefore, R^ a (x,y) is a kernel of a Hilbert-Schmidt operator for any fixed 
A 6 (—00, —1] and a > |. 

It remains to show that K^ a (x,y) is uniformly bounded in the limit A — > —00 (k —> —00) 
for any i,|/6K. Note that 

u (z.A^m (y,fc)j -[ u+{x . k)2U - {y . kh u+fakhu-fak)! 

and a similar formula for u~ (x; fc)[<7iU + (?/; k)] . By Proposition [5j this matrix grows linearly in 
k as k — > —00 for any x,y £ R. On the other hand, fea+(/c) grows at least linearly as \k\ —> 00, 
which implies the A-uniform bound ([66|h 
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To prove bound (|67p. we can see from the linear growth of u^(x;0) as x — > =Foo that 
K^ a (x,y) is a kernel of a bounded operator from L 1 (M) to L°°(R) for any a > 1. The 
mapping is also bounded as k — > — oo. □ 

Let P a .c.{Ti) '■ L 2 (R) i — y L 2 (M) be the orthogonal projection operator to the continuous 
spectrum of %. We recall the Cauchy formula, 

(68) e- im P a . c .(H)f = ^- (| 1 + ^°°) e- ux [R+(\) - it^(A)] fdX, 

where the integral is understood in the norm of the mapping from L 2 a (R) to L 2 _ a (R) for a > 3 
The interval (— oo, —1] for A can be parameterized by (— oo,0] for k using the substitution 

kdk 



2 1 



A = -Vl + A: 2 => dX 



These representations are used for the derivation of linear dispersive decay estimates for the 
semi-group e~**^P ax .('%). 

7. Linear estimates for the operator % 
We shall need two preliminary results, which will be useful in our arguments for this section. 

7.1. Preliminaries. The first result that we need is the Christ-Kiselev lemma. We actually 
state a version due to Smith k, Sogge 



Lemma 1. Let X, Y be Banach spaces and K, : L P (R; X) — > L q (M., Y) be a linear operator such 
that tCf{t) = f_ K(t, s)f(s)ds. Then, the operator 



(69) £/(*)= / K(t,s)f(s)ds, 

Jo 

is bounded from L P (M;X) to L q (R, Y), provided p < q. Moreover, there is C PtQ > such that 



\\fc\\LP(M;X)-+L<l(m.,Y) < Cp,<2 \\JC \\ Lp(M.;X)->L9 (M,Y) • 

The second lemma is a technical statement, which is complementary to Lemma HJ when the 
condition p < q is violated (most notably when p = q). This is stated for the Schrodinger 
operator — d 2 + V(x) by Mizumachi (Lemma 11 in [IT]), but it applies equally well to an 
arbitrary self-adjoint operator C. 

Lemma 2. Let £ be a self-adjoint operator and P a , c (C) be a projection to the absolute contin- 
uous spectrum of C. Let g{t,x) = g\{t)g2{x) and define the function 

■ l-OO 

(70) U(t,x) = -= / e- rfA #i(A) ([Rt(X) + R C (X)] g 2 ) (x)dX. 

Then, we have 
U(t,- 



rt / r0 r°°\ 

2 / er^ c P a . c XC)g(s,-)ds+( / - / ) e~^ c P a . c .(C)g(s, -)ds. 

JO \J-oo Jo / 



We use the resolvent analysis of the Dirac operator T~L to derive some linear estimates, which 
are used in the proof of the main theorem. 
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7.2. Mizumachi estimates. We refer to Mizumachi's work [T7] in the context of the one- 
dimensional NLS equation, which was used in our work [12] in the context of the discrete 
NLS equation. These estimates are developed to control quadratic nonlinearities in the time- 
evolution equation (|41|) . which have fast spatial decay. Thus, the challenge here is to achieve 
L\ temporal decay, in the presence of the exponential spatial decay. 

Lemma 3. Fix a > |. There is C a > such that 



(71) 
and 

(72) 



- a e- m P a .c.CH)f\\ L ^ L? <C a \\f\\Ll 



(x)- a r e -^-^P a . c .(^)F(r,.)dr 
J o 



<C a ||(x)«F|| L i L? . 



Proof. The proof of Lemma [3] proceeds via analysis of the contribution of the high energy part 
and the low energy part. 

Let x( x ) be an even C°° function with x( x ) = 1 f° r \ x \ < 1 an d x( x ) = for \x\ > 2. Fix 
M > 2, let xm( x ) = x( x /M) and decompose 



Pa.c.mf = XMe~ im P a . c m)f + (1 " XM)e- UH PaAK)f- 



an 



In order to show (|7T|) . we need the following two estimates 

(73) ||(1 - XM)e- m Pa,Xn)f\\ L ~ L t < C\\f\\ Ll 

(74) || < x >-* XMe- m P a AU)f\\ LTLl < C\\f\\ Ll 

Combining bounds ([73]) and ([71]) . we complete the proof of estimate ([7T]) . Bounds ([73]) and 
(|74p are proven in the following two subsections. 

The proof of estimate (|72p is based upon Proposition [8] and Lemma [2] By Lemma [2] we can 
write (with C = %) 

J* e-^ n P a . c .{U)F{T, -)dr = hj + \ (T ~ J_) e-^ n P a . c .m)F(r, -)dr, 

where U is the function defined by (|70p . Let us first control the last two terms. Since they are 
similar, we only need to control one of the terms. By the estimate (|7ip . we have 



(x)- a e- itn 



/>oo 

/ e tTn P a . c .(n)F(r,.)dr 
Jo 



< C 



/ e lTn P a . c .(H)F(T,-)dT 
Jo 



< C\\{x) a F\\ LlLh 



where in the last step, we have used the dual estimate to (I7ip . In order to control the U term, 
we observe that the set of all functions {g\{t)g2{x) : g\ G P\-,g2 £ L^.} is dense in L\.L\. The 
estimate that we need follows from 



<*)- 



e- itx g\(X) [R+(X) + R^(X)] g 2 dX 



<C\\ gi \\ L 2\\(x) a g 2 \ 



The left-hand side is controlled by Minkowski's inequality and Plancherel's theorem in the time 
variable, 

\\(x)- a \\g\(X) [R+(X)+R n (X)] salltfJU- < C||gi(A)|| L 2 sup ||i2±(A)|| L ^ L 2_ J(x) a g 2 \\ L ~. 



Using bound ([67]) of Proposition[8]for any a > 1, we bound the last expression by C\\g\ \\ L 2 ||(x)<72||li , 
which completes the proof of estimate (I72p . □ 
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7.2.1. Proof of (|73p . Using the Cauchy formula (|68p for 

g x>t (X) := (l-XM(\))e- m P a .c.(n)f, 

we can see that for each fixed value of x, this function is a multiple of the Fourier transform of 
the function 

g x (X) := (1 - xm(\)) {[B+(X) - R n (X)] f) (x), 
evaluated at t. Therefore, by Plancherel's theorem, we have 

\\(l-XM(X))e- itH P a . c .(H)f\\ Ll = C\\g x \\ Ll . 

It is sufficient to control 



sup 11(1 - XM(X))R±(X)f(x)\\ L] { <C\\f\\ L z, 



which we will do next. By iterating the resolvent identities, 

R h = Rq- RhVRq = R - RoVRh, 
we get the representation formula 

(75) Rh = Ro~ RhVRq = Ro~ RqVRq + R Q VR n VR . 



where Rq is the resolvent of the free Dirac operator D defined by (jlip . For the first term, we 
have 

sup ||(1 -XM(X))R±(X)f(x)\\ L 2 
By symmetry, it suffices to consider only positive values of A, whence we need to control 

\R^(X)f(x)\ 2 dX. 



sup 

xGRJM 



We compute the resolvent Rq(/j.) 

R Q ( P ) = (D-ny i = (i-d 2 x -n 2 )- 1 



-id x +n -1 
-1 id x + n 



for fi ^ °~(D) = [— oo, —1] U [l,oo]. By analytic continuation, we may define the resolvent even 
for values on the spectrum of c(D). Since we need such a formula for values of fi £ (M, oo), it is 
convenient to introduce a change of variables /x = \fk 2 + 1. Note that dfj, = k(k 2 + l)~ l l 2 dk ~ 
dk and the interval of integration becomes (V M 2 — 1, oo). Now, since the resolvent operator 

(—d 2 — k 2 ±0) _1 is given by a convolution with the explicit kernel 6 1 , it is clear that R^(fj,)f 
is a linear combination of convolution operators with kernels 

e ±ik\x\ e ±ik\x\^/f,2 _|_ ]_ 



(76) 



a ±ifc[x| 



sgn(x) 



k ' k 

We shall consider the first type of operators, the second one has a stronger decay, while 
the third one is basically the same as the first one. By Plancherel's theorem applied to the 
functions f(y)Xy<x and f(y)Xy>x, we have 

2 



VM' 2 -1 



< 2 



VM 2 -1 



f{y)dy 



2 




+ 


/ 

J X 



dk 



±iky 



f(y)dy 



dk<C\\f\\%. 
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Similarly, we estimate the contribution of the second term RqVRq in the expansion ()75|) . 
Again, we have to deal with different terms of the convolution operators, but the hardest one 
is again e ik ^sgn(x). We get 



pOO poo f 

/ / e ±lk \ x - y hgn(x - y)V(y) / e^^sgn^ - z)f{z)dzdy 

JVM 2 -1 J-oo J 



dk 



< C\\V\\ 2 L1 sup 



VAP-l 



dk < c\\v\\l\\f\\l l: 



where in the first inequality, we have applied Minkowski's and at the second inequality, we have 
applied our previous estimate. 

In order to estimate the last term in (|75p . we use bound (|66p of Proposition [8] and get 



VM 2 -1 



±ik\x-y\ 



sgn(x - y)V(y)[E±{y/l + k?)V{R±{sfl + &)f)(y)dy 



dk 



poo 

< C\\ < x > a Vf L2 / < y >- a R±(Vl + k*)VR±(Vl + k*)f 

x JsJM 2 -\ 

roo 



dk 



poo 

<C\\<x > a Vf L2 sup / \R±(^l + k 2 )f{y)\ 2 dk 

<C\\<x><*V\\%\\f\\%. 
This concludes the proof of (fT3j) . 



7.2.2. Proof o/ (HU). We shall prove that 

(77) sup < x >- 3 / 2 [| XM (A)(i^(A)/)(*)|U < C\\f\\ Ll , 



which implies bound (|74p by Plancherels' theorem and Cauchy's formula (|68p. To prove (|77p 
for A < — 1, we use representation (|62p and write explicitly 



XM(\)(Rt(\)f)(x)\\ 2 L2 



-1 



X 2 M (X)\(Rt(X)f)(x)\Ux 



\i(x,k)\ 2 \k\dk 
-Vm*=1 ^k 2 \a+{k)\ 2 y/l + ^' 



where 



f(x,k):=u+(x;k) I [a 1 M-{y-k)\ 1 f{y)dy + vT(x-k) I [a^ (y; k)] T f(y)dy. 



For definiteness, let us assume that x > 0. We represent 



u-(y;k)} T f(y)dy = f '[am' (y;k)] T f(y)dy + f [a_, l]f{y)e- ik »dy 

JO J-oo 

+ f {[a im r{y;k)] T -[^-A]) f(y)e- tky dy^h + I 2 +h 

J — oo 



and 



[a 1 u + (y;k)] T f(y)dy 



POO 

[a+, l]f{yy ky dy + / {[axm+iy; k)f - [a+, 1]) /(y)e^dy 

J X 



h + h- 
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Using Proposition [5] and Cauchy-Schwarz inequality, we have 

< llu-Oj^ll^^ll/ll^^c^) 3 / 2 !!/^, 
|j 3 | < iim-^fej-^a-fn^^ii/ii^^qK^ViUyii/ii^, 

\h\ < ||m+(-;A ; )-[l,a + ] T || ii(R+) ||/|| L , < C\\(x) 3 V\\ L ~\\f\\ Ll . 
The estimates for I3 and I5 follow from the bound 

(78) ||m+(.;/c) - [1, a+] T || L 2 (R+) + Urn" (•;*:) - [1, a_] T || L 2 (ffi _) < C\\(x) 3 V\\ LS > 

which we prove now. We need only control the first term, the other one is controlled in a 
similar matter. 

By the formula (|48|) . for all x £ R and all £ R near A; = 0, there is C > such that 

|G+(x;A;)| < C(x>. 

By Proposition [5j for all x > 0, there is C > such that 

\m + {x,k)\ = |u + (x,A:)| < C. 
Thus, by the integral equation (f47l) . we get for all x > 0, 



OG 



|m + (.r:A-)-[l.n + ]'| < Cj (x - y)\V{y)dy < C\\(rfY\\ L> / (:)— -^/.: 

< C||(x) 3 y|| L? c <x>- 1 



(a; + z)' c 



This computation completes the proof of the first inequality in ([78]) . 
On the other hand, for any finite M > 1, PlanchereFs theorem gives 

(l^| 2 + |/4| 2 )^<C||/||| 2 . 

v/M 2 -l 

Since ka + {k) is bounded away from zero as k — > and |a + (/c)| > 1, we obtain 

which concludes the proof of bound ()77p and hence of bound (|74p . 
7.3. Strichartz estimates. We use the following standard definition. 

Definition 1. VFe say that a pair (q,r) is Strichartz admissible for the nonlinear Dirac equa- 
tions if 

2 11 

q > 2, r > 2 and - + -<-. 

q r 2 

In particular, (q,r) = (4, 00) and (q,r) = (00, 2) are end-point Strichartz pairs. 

Lemma 4. Let (q,r) be a Strichartz admissible pair, s > 0, and e > 0. Then, there are 
constants C £ > and C > such that 

(79) l|e-^Pa. c .(WIL4£~ < C £ \\f\\ H s /i+s , 

(80) ||e- im P a . c .(^)/|U ri?l < C\\f\\ m , 



(81) 



f e -i(t-T)Hp ac {H)F ( Tr)dT 
Jo 



< C\\F\ 

-'t J - 1 x ] ]J ^t 



L?° HlnLf Lt 
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Proof. Let us first comment on the estimates (|80p and (|8ip . It is easy to see by the self- 
adjointness of H that (|80j) is trivial for s = 0. We easily extend to all integer values of s by the 
observation that d x behaves like H and commuting 7~L with e~ lt ^P a , c , (H). This is made precise 
in formula (|86p below. We then conclude by interpolation to obtain (|80p for all nonnegative 
values of s. Regarding (|8ip . it follows by an easy application of Lemma Q] combined with the 
dual estimate of (fT9j) . 

Thus, it remains to show (|79p . We will in fact deduce this Strichartz estimate for the 
perturbed Dirac operator ~H by using the corresponding result for the free Dirac operator D, 
in addition to the weighted estimates in Lemma [3j This is in essence the approach taken by 
Rodnianski and Schlag, [22]. Let us first record the Strichartz estimates for the Dirac operator 
D 

(82) We^fh ^ ' ' ' 



< c s II f\\ H s(*,r), s(q,r) 



1 1 



for all Strichartz admissible pairs (q,r), so that q > 4 + 5. This of course looks exactly the 
same as the estimates that one gets from interpolating between (|79p and (|80p . We refer the 
reader to recent work of Nakamura-Ozawa, [19] (more specifically Lemma 2.1 with 9 = 1, A = 
3/2, n = 1) for a reference for this result. Note that this result would not extend to the full 
range q = 4, r = oo, unless we are willing to replace the L°° by the Besov space 2 (which 
we are avoiding for the purpose of simplicity). In order to extend this to the useful endpoint 
q = 4, r = oo, we must introduce slight loss of smoothness, so we have 



(83) 



-itD 



/II 



<a 



Fix now e > and take a test function / 
D + V(x), we may write 



P a .cXU)f G # 3 / 4+e . Recall that since U 



-itU 



f = e - itD f - i f e - i{t - s)D Ve- isn fds. 
Jo 



Furthermore, we may write the symmetric matrix V(x) as the product of V\(x) and V^x), where 
both V\{x) and V^(x) are C^-smooth and have fast decay at spatial infinity. For instance, one 
may pick V\{x) = V(x)(x} 10 and V^x) = (x)~ 10 Id. We have 



e- m Pa.c.CH)f\\ L f L ~ < \\e- ltu f\\ LfL ~ + 



itD 



< a 



H 



3/4+e + 



- i ^ D V 1 V 2 e- isH P a . c .{%)fds 
e-< t -^ D V 1 V 2 e- isH P a . c .(n)fds 







At this stage, in order to estimate the second term, we will use Lemma [TJ Let K(t, s) = 
e -i(t~s)Dy^ considered as acting between LlH^ A+£ to L\U£ . The Duhamel's term that we 
need to estimate is 



M(t) 



K(t, S )V 2 e- lsn P a . c .(n)fds = KV 2 e- UH P a . c .(H)f 



an 



where K is defined by (|69p . It follows from Lemma Q] (since q = 4 > 2 = p, this lemma can be 
applied) that 



\M 



\l*L°° < C \\ K \\^ t H^ i+ ^LiL^ V2e tmP a.c(^)f\\ L 2 H i 



3/4+ £ 
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We need estimate then the operator norm \\K\\ 2 „3/4+ e r4roo - We have by ([83 

/OO 
-00 

We will show that for < s < 1, 



e ifl£ ViG(s, -)ds 



H 



3/4+e 



(84) 



e* si ViG(s, -)ds 



< C s .v, I Gil 



and 

(85) \\V 2 e- im P a . c .(H)f\\ Ltm < C V2 \\f\\m, 

which implies what is needed. Indeed, for s = 3/4 + e, we deduce 



IMI 



L | L ~ < C Vl \\V 2 e- im P a . c XH)f\\ L2HS/ 4 + e < C VuV2 \\f\\ H 



3/4+e • 



It thus suffices to establish (|84p and (|85p . By interpolation, it suffices to check both only for 
s = and s = 1. The statements for s = in fact follow from the corresponding arguments for 
s = 1, so we concentrate on s = 1. For (|84p . (observe that d x e ltD = e ttD d x ), we have 



e isZ ViG(s, -)ds 



< 



e isZ ViG(s, -)ds 



+ 



LI 



LI 



By the dual estimate to (I82p (recall s(oo,2) = 0), the right-hand side of the last inequality is 
estimated by 

c(\\ViG(s,-)\\ L}Ll + mvMs.-MLiLi) < emu- + \\vI\\l~)\\g\\ l}h1 . 



This is the proof of ([84]) . 

Next, we need to deal with derivatives in the estimates for the perturbed evolution. From 
the formula D = % — V(x), we have the equivalence 















C) 






+ 


C) 




H 1 






L 2 





L 2 

which will be used repeatedly in the arguments to follow. Regarding (|85|) for s = 1, we use 
(1861) to obtain 



,2 . 



ll^[F 2 e-^/]|| L?L , < WVZe-WfWqq + \\V 2 e- im Uf\\ L , Ll + \\V 2 e- m f\\ L , L 

Now, since |V^'(x)| + |T^(x)| < (x)~ 10 , we estimate the last three quantities by 

C\\(x)-H Ll ^)- 5 e~ UH f\\ L ^ L 2 + \\(x)- 5 e~ un Hf\\ L ^ L 2 + \\(x)- 5 e- m f\\ L ^ L : 

< leqC(\\f\\t> + \WWz3) < C\\f\\ H . 



where bound ([TTj) and Holder's inequality are used. This computation establishes ([85[) and 
hence Lemma 0[ □ 
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7.4. Additional estimates. Mizumachi estimates and Strichartz estimates admit a number 
of useful corollaries. 

Corollary 1. Let (q, r) be and admissible Strichartz pair such that q > 4 + 5. For each 5 > 0, 
there is Cs > such that 



(87) 



1 1 e iTH P a . c {U)F(r,-)dr 
Jo 



<c q 4FL„, 



Li 



. . 1 1 1 

S{q > r) = 2 + q-r-> 



where (q',r') are duals of(q,r). 

Proof. The result is obtained from ([79]) and (]80p by duality and interpolation. 
Corollary 2. Fix a > 2. There is C a > such that 



□ 



/■* e -i(t-r)«p ac( ^ )F(T) . )dr 
JO 



<C a (\\{xTF\\ LlLl + \\{xrd x F\\ LlL ,). 



Proof. Due to the density of {<7i(i)<72(a0 : 51 6 Lu92 6 -^i} in L\.L\, it will suffice to show 



(89) / e^-^P a . c .(^)<?i(r)52dT 
Jo 

By Lemma [H we need to show that 

/oo 
e -«t-T)Hp ac{n)gi ( T)g2dT 
-oo 

By d79|) and (HO]), we have 



LfHinLfL? 



<C|| 5l || z? (||(x) Q 52 || z i + ||(x)^ 2 || z i) 



LfHinLfLS> 



<C\\g 1 \\ L? (\\{xrg 2 \\ Ll + \\(x) a d x g2\\ Ll ) 



e lTn P a .c.(H)gi(T)g2dT 



< 



LfHinL*LS> 



e lTn Pa.c.m)gi(T)g 2 dT 



Hi 



Again, one may convert one derivative to T~L — V by the equivalence (J86J), whence we further 
estimate by the dual of ([71]) . 

COO 

e lTn P a .c.(H) gi (T)g 2 dT 

> #1 



< c 



e lTn P a .c.(n)gi(T)Hg 2 (-)dT 



+ 




LI 





V I e^ H P a . c X%) gi {r)g 2 {-)dr 



LI 



<C\\gi\\ Ll (\\(x) a Hg 2 \\ Ll + \\(x) a g 2 \\ L i) 
<C\\gi\\ L 2(\\(x) a g 2 \\ Ll + \\{x) a d x g 2 \\ Ll ), 
which is the desired estimate. 

Corollary 3. Fix a > 2. There is C a > such that 



□ 



(90) 



<*)- 



-i{t-r)H 



P a .cXn)F{r,-)dr 



< C\\F\ 



More generally, let (q, r) be an admissible Strichartz pair. Then, 



(91) 



(xY 



f e-^-^P a _ c XH)F{r,-)dr 
Jo 



< C n \\F\ 



L q t 'wy> 



where (q',r') is a dual pair. 
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Proof. The proof of (|90p is by averaging the estimate (|7ip . More precisely, using the triangle 
inequality and estimate (fTTI) yields 



(x)- a / e-^-^Pa.c.WFir,-)^ 



< C 



/oo 
\\(xr a e- i(t - T)H P a ,.(n)F(T,-)\\ LT , Li dT 
-oo 



-oo 

"OO 



/oo 
\\F{T,-)\\ Ll dT = C\\F\\ LlLl . 
-oo 

For the proof of (I9ip . we use Lemma [TJ It will suffice to bound the operator 



TF(t) = (s ) 



for any fixed xq G K. We have, by ([7T]) 



|TF|| L? < 



a-iffl 



-i{t-r)H 



e iTn F{r, -)dr 



F(t, -)dr 



< C 



e irW F(r, -)dr 



By Corollary [[J we bound the last expression by 



as stated in ([9T]) . In the last step, we have used that if 4 < g < oo and 5 <C 1, then ^+<5 < 1. □ 

Corollary 4. Fix a > 2. There is C a > such that 

(92) IK^^^e-^Pa.cCW)/!!^ < C a ||/||m, 



(93) 



anc 



(94) 



Jo 



<C\\F\\ LlHi , 



(x)- a d x f e-^- T ^P a AU)F{T,.)dT 
Jo 



<C(\\(x) a F\\ L ^ L? + \\(x) a d x F\\ LrL 



■ 4' 



Proof. The proof of the estimate (|92jl is based again on the equivalence (J86J). Since Ti commutes 
with all functions of T~L (by the functional calculus), we have from (|86p and (|7ip 



(x)-«d x e-™P a . c .(H)f\\ LTq < \\(x)- a e -^P a . c .(n)nf\\ L?L 2 + \\(xr a Ve-^P a . c .f\\ LTL2t 

< C(\\nf\\ L2 + \\V\\ L ^\\f\\ L ,)<C v \\f\\ m . 



The proof of the estimate (|93|) is by averaging. Indeed, by (|92j) . 

(x}~ a d x I e-^-^P a . c .(^)F(r,-)dT 
J o 



< C 

< C 



/oo 
||(x)- Q a :E e-^-^P a . c .(^)F(r,-)|| LS o L? cir 
-oo 



— oo 
oo 



F(T,-)\\ m dr = C\\F\ 
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For the proof of the estimate ([94"]) , we apply again the equivalence ([86]) and then we use ([72]) . 

We have 

-t 



(x)~ a d x [ e- i ^ H P a . c X'H)F(r,-)dT 
Jo 

c)- a f e- i ^ u P a . c .{U)'HF{T, -)dr 
Jo 



+ 



<x>- 



W{x) f 
Jo 



-i{t-r)H 



P a .c.(U)F(T,-)dT 



<C(\\(x) a HF\\ LlL2t + \\(x) a F\\ Llxq ) 
<C(\\(x) a F\\ LlL , + \\(x) a d x F\\ LlL ,). 
This concludes the proof of the corollary. 

8. Proof of the main theorem 
We first formulate the solution and the nonlinearity spaces. Let 



□ 



•2 + 



+ \\Y\\ L?Hl , \\Y\\ X2 :=\\{x)- a Y\\^ 
and ||Y||x := HYH^ + ||Y||x 2 - The nonlinearity space is defined via the norm 

inf IIF, + (\\{xr¥ 2 \\ LhL 2 + \\(x) a d x ¥ 2 \\ LlL ,). 



\F\\m := inf ||Fi||ri Hl + 

F=Fi+F 2 "^t^x 



Consider the Cauchy problem for the inhomogeneous linear equation, projected along the 
absolutely continuous spectrum of H 



(95) 



if = OT + P,,(^)F, 
Y(0)=Y = P a . c .("H)Y 



When one interprets correctly the results of the dispersive decay estimates (Section 7) in the 
notations above, we get that a solution to the Cauchy problem (|95[) satisfies 

(96) \\Y\\ x <C(\\Y \\ m + \\F\W). 

For the proof of the main theorem, we need to show the existence of small solutions for the 
system of two (scalar) ordinary differential equations f)42[) for u and 9 coupled with the partial 
differential equation (j44p for Y. 

Since the right-hand side of equation ([44]) is not projected to the continuous spectrum of H, 
we decompose 

(97) Y = au + Z, a = (u ,Y) L2 , (u ,Z) L2 =0, 

where u is the eigenfunction of % for eigenvalue ujq. Substituting ([97]) into ([4"I]) . we obtain 
the system of equations 

ia = u a + (u , e~ %e F) L 2 , 
iZ = HZ + P a . c .(H)e- ie F. 

We now set up our problem as an iteration scheme, where we look for a fixed point in a small 
ball in a normed space. More precisely, this space is composed of all quadruples (u,6,a,Z), 
equipped with the norm 



(98) 



M\l\ 



+ \\ a \\L 2 t nLf 



+ IIZI 



X ■ 



\\(co,9,a,Z)\\ z :-- 

Note that the elements of the corresponding set are subject to the appropriate initial conditions 
w(0) G (ojq, ojq + e), 9(0) = 0, a(0) = (u , Y(0)) L 2, Z(0) = P o . c .(ft)Y(0). 
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First, observe that the matrix in front of the variables Co and 9 — u in (|42p has the form 
(99) 

Due to the smallness of 



Re(a w U,U-Ui} L2 Im(^U,Ui) L2 
Im(^U,Ui) L2 Re(a w U,U + Ui) L2 



i^liuiii.w + odiuxii^ 



l|Ul||L t °°L2 < \W\\Lf + ||Z|| L oo L 2 

and the non-degeneracy condition (|35p . we may conclude that the matrix (|99p is invertible. 



(Note that ||Z|| £00^,2 is a part of the norm ||Z||x, which is kept small in our fixed point argu- 
ments.) 

Next, we show that the quantities H^H^i and — w[|x,i are under control. Indeed, due to the 
invertibility of the matrix, the form of (|42p . and the quadratic nature of F\,F2 (Proposition 
HJ), we have 



W rl + 



| L i < CdlFilLi + 11^11^) <C 
CII < a; > 2a U 



-.1 < a; >" 



|Y(x,t)| 2 |U(x)|dxdt 



l LsoL 2<C(||a|| 2 ? + ||Z| 



It follows from this bound that 
(100) 



oj - w(0)||loo + ||0 - f uj(s)ds\\ L oc < C (\\a\\ 2 L 2 + ||Z| 
Jo \ t 



Since uj € Lj and — w(0)||£,« is small, there exists eo > e and := lim^oo cj(i) such 
that Woo G (wo,o;o + eo) if w(0) G (wq, wo + e). Similarly there exists 6*00 G R such that 



lim (e(t) - [ 
V Jo 



uj(s)ds 



Now, we control the quantity ||a|| i 2 ni oo. It follows from the symplectic orthogonality condi- 
tions (EHD that 



(u ,Ui) L 2 = Re(u 



Ui) L 2 +ilm(u - ,Ui) L 2. 



JU|| L 2'-'^ ' — * u ||^U 
By Proposition [21 for any a > 0, there is C a > such that 



< x > a (u 



U 



U 



+ 



Li 



< x > (u 



-) 



< C\uj — UJq\. 



LI 



Therefore, we obtain 

\H\lI = \\( u o, y )l%\\l1 = IK u 0,Ui) L 2|| i2 

< C||w-wo||£»|| < x >~ a Ui|| ig0£ ,2 < C(e+ ||w -w(0)||£|o)[| < x >~ a Y\\ L ^, L 2, 

where e+ \ \uj — (jj(0)\\l5° is small due to smallness of e and the bound (jlOOp . Similarly, we obtain 

IMUf = IK U 0, Y) L 2 || L oo = ||(u ,Ul) L 2|| L oo 

< C\\uj -wo||l°°||Ui|| l? o L 2 < C(e+ 11^-^(0)11^00)1^11^00^2. 

Finally, it remains to estimate the quantity ||Z||x- Due to our construction, we have Z = 
^o.c.(^)Y, so that we may apply the linear estimates ([96]) . The nonlinearity P a . c .(H)e~ ie F in 
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the residual equation (|98p has two parts. The first part satisfies 

||Pa.a(ft)e-*(«W? w U + (9 - u>)U)\\ LlH x < C(\\u\\ Ll + \\6 - a;|| x i)(||U||^ + H^UH^) 

< c(\\a\\l* + \\Zf x ). 
Roughly speaking, the second (nonlinear) part 

G := P a .c.W)e- ie (N(U + Ye ie ) - N(U)) , 

consists of the two terms 

G~ YU 2p + Y 2p+1 , 
where Y is controlled in the X-norm by 

\\Y\\ x <C(\\a\\ qnLr + \\Z\\ x ). 

Note that 

|G(a:,t)| + \d x G(x,t)\ < C(\Y\ + |^Y|)(|U| 2 f + \d x XJ\ 2p ) + C(|Y| + |^Y|)|Y| 2 f , 
We need to control the quantity 1 1 C | |_a/" in terms of ||Y||x- We have 

\\G\W < C\\<x >« (|Y| + |^Y|)(|U| 2 p + \d x XJ\ 2p )\\LlL? + <7||(|Y| + I^DIY^H^ 
< C (\\(xr a Y\\ L ^ + \\(x)- a d x Y\\ L ^ L i) || < x > 2a (|U| 2 ^ + \d x U\^)\\ LlLr 
+C||Y|| Lt » H i||Y|| 2j L . 

It is now easy to close the argument in the norm ||Y||x- Indeed, by Sobolev embedding for 
any e > 

||Y||x-l« < C\\Y\\ L ^ Hl ,^ < C\\Y\\ X . 

We also have ||Y||£4 L oo < ||Y||x (by the definition of || • ||x) and hence, for p > 2, by the log 
convexity of the L q norms, we have 

\\ Y \\ L f L ~ < II y IIl4l ? II y IIl-l°° <C|I Y IU- 
All in all, combining the estimates for ||G||_v with the estimates for ||Y|L2 Proo , we obtain 

t X 

\\G\W <C\\<x > 2a (|U| 2 p + \d x U\ 2 n\\ L iL r \\nx + C||Y||^ +1 . 

By Proposition O there is C > such that 

|| < x > 2a (|U| 2p + |^U| 2 f)|| L i Lr < C\\u - wo|U- < C(e + \\u - w(0)|U«). 

Since the last term is small due to the smallness of e and the bound (jlOOp . the fixed point 
argument is closed, and the proof of Theorem [1] is complete. 
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